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ABSTRACT 

We present calculations on the formation of massive black holes with 10 5 M 0 at z > 6 
that can be the seeds of supermassive black holes at z > 6. Under the assumption of 
compact star cluster formation in merging galaxies, star clusters in haloes of 10 s ~ 
10° M 0 can undergo rapid core-collapse leading to the formation of very massive stars 
(VMSs) with ~ 1000 M 0 which directly collapse into black holes with similar masses. 
Star clusters in halos of > 10 9 M 0 experience type-II supernovae before the formation 
of VMSs due to long core-collapse time scales. We also model the subsequent growth 
of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently 
re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small 
radii is the main driver for accretion of stars onto black holes. As a result, more than 
ninety percent of stars in the initial cluster are swallowed by the central black holes 
before z = 6. Using dark matter merger trees we derive black hole mass functions at 
z = 6 — 20. The mass function ranges from 10 3 to 10 5 M 0 at z < 15. Major merging 
of galaxies of > 4 x 10 s M 0 at z ~ 20 successfully leads to the formation of > 10 5 M 0 
BHs by z > 10 which can be the potential seeds of supermassive black holes seen 
today. 

Key words: stars: black holes - quasars: supermassive black holes - galaxies: star 
clusters: general - galaxies: formation galaxies: high-redshift 


1 INTRODUCTION 


Recent observations of high-redshift QSOs at 2 > 6 sug¬ 
gest supermassive black holes (SMBH) of Mbh U 10 9 Mq 
form on short time scale < 1 Gyr (e.g ., |Fan et alj 12001 . 


200f :lMortlock et al. 11201 iHKotrnerid v fc Holhoi til: Wu et alj 


20151 ). The most distant QSOs are observed at 2 = 7.1 which 
corresponding to an age of the Universe of ~ 800 Myr. The 
formation mechanism of SMBHs on such short time scales 
has not been understood yet. Cosmological simulations show 
that B Hs can grow up to ~ 1 0 9 Mq via Bondi-H oyle ac¬ 
cretion dDi Matteo et al.ll2008l . I 20 I 2 I : iLi et alJl2007h . While 


there are uncertainties associated with the modelling of the 
accre tion efficiency at the sub-grid level dBooth fe Schava 
l2009h . it appears that the main crucial assumption is that 
massive black holes (MBHs) of 10 5 M@ are seeded in haloes 
of ~ 10 10 Mq within the simulations. These simulations in¬ 
dicate that once MBHs form in the first galaxies, they can 
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grow to SMBH via gas accretion at Eddington rate by z ~ 6. 
The main obstacle is thus referred to forming MBHs. 


One natural scenario for the formation of massive BHs 
is the Eddington-limited growth of a few 100 Mq stellar 
mass black hole via gas accretion that is the remnant of a 
Population III star. However, due to stellar feedback, the gas 
accretion rate is s ignificantly suppressed and far below the 
Eddin gt on-limit d Johnson et all |2007; Milosavlievic e t al.l 
l2009bl lal: I Alvarez et aI.ll2009l ; |Park fe Ricottill201ll . l2012l'l . In 
addition, cosmological simulations show that Population III 
star remnants mostly resided in low-density environments, 
hence the gas accretion rate is much smaller than the Ed¬ 
dingt on rate even without stellar feedback dAlvarez et aid 

l20oJl . 


Another proposed mecha nism is the direct collapse 
of su per massive stars (SMS) ( ReesHl978l : iBromm fc Loebl 
2003). If the formation of hydrogen molecules in the 
first galaxies is suppressed by external UV radiation 
in the Lyman-Werner bands, gas clouds cannot frag¬ 
ment. As a result, the clouds cannot form low mass 
stars but form SMSs of ~ 10 5 Mq due to high-gas 
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2 Yajima et al. 


accretion rates of ~ 0.1 — 1.0 M 

7) yr 1 (i Omukai et al. 

20051: Lodato & Nataraianl 200(1. 

20071: Begelman et al. 

200fi!: Soaans & Silk 200fiJ: Inayoshi & Omukai 

2012i: 

Agarwal et al.l 20121. 120131: Hosokawa et al. 2013 

). This 


scenario is also supj 

sorted by detailed numerical simulations 

(Regan & Haehnelt 

2009: Regan et al. 

20141: Maver et al.1 

201CJ; 

Choi et al. 2013: Latif et al. 20131; 

Johnson et al. 20141: 

Inavoshi et al. 20141: Agarwal et al. (2014|). One of the main 


requirements of this scenario is a strong UV background ra¬ 


diation with J 21 > 100 in units of 10 21 erg s 1 cm 2 


Hz 


and very low metallicity (Z < 10 ~ 5 Z Q ) (Omukai ct al. 


2005|j200^jDjjkstraetaljlgOO^A garwal feKhochfarl 20151 : 


Sueimura et al. 2014 : Inavoshi fe Tanaka 2014h . Note that, 
however, the critical value of J 21 is still under the debate. 
It sensitively d epends on the environments via, e.g., the 
shap e of SEP ([ Sugimu ra et al.l 120141 : lAgarwal fc Khochfad 
120151 : lAgarwal et alJl2015hh . 

Recently iMaver et al.l (l2010l l showed in high-resolution 
simulations that another suggested path to the formation of 
MBHs via major mergers of galaxies could work. Their sim¬ 
ulations showed that a large amount of gas in disc galaxies 
falls down to the galactic centres due to angular momen¬ 
tum loss by tidal force, and massive high-density gas clumps 
form. They suggest that the gas clouds directly collapse into 
MBHs. 


In addition, llnavoshi et al.l d2015h suggested that SMSs 
could form in merging primordial haloes due to collisional 
dissociation of H 2 in the compressed high-density regions. 
The outcome of the collapse depends sensitively on the equa¬ 
tion of state of the gas. If radiative cooling quickly occurs, 
a collapsing gas clo ud fragments and forms a dense star 
clust ers instead (e.g., [Regan fc HaehnehJl2009l : lFterrara et al.l 
1201311 . 


Due to the gravothermal collapse, such dense star clus¬ 
ters can cause the merging of massive stars and runaway 
growth, resulting in the fo rmation of very m assive stars 
(VMSs) of ~ 1000 Mq fe.g.. iFre itag et akll^OOfSi h Direct nu¬ 
merical simulations bv IPortegies Zwart fc McMillanl d2002l l 
support this view and show that compact dense star clus¬ 
ters cause core-collapse and make VMSs with i n the typical 
lifetim e of massive stars (see also, iFuiii fc Portegies Zwartj 
I 2 OI 4 I 4 . VMSs can result in BHs with almost the same mass 
ijHeger et al]l200H h Subsequently the BHs can grow via ac ¬ 
cretion of stars with ti d al dis ruption events fe.g. lReeslll978ll . 
iDevecchi fc Volonteril (2009h analytically modeled the for¬ 
mation of VMSs in first galaxies via disc instabilities, and 
showed a fracti on ~ 0.05 of first gala xies at z ~ 10 — 20 form 
VMSs (see also IDevecchi et alJl2012h . 

Very recently iKatz et al.l (120151 1 followed the merging 
processes of metal enriched haloes in detailed cosmologi¬ 
cal hydrodynamics simulations and showed the formation of 
high-density gas clumps which are potential sites of dense 
star clusters leading to the formation of VMSs. In this work, 
we investigate the possibility of the formation of MBHs via 
the formation of VMSs and subsequent growth by stellar 
relaxation processes in merging galaxies at 2 > 6. 


Our paper is organized as follows. We describe our 
model in Section [2] In Section [3j we show the masses of 
VMSs as a function of halo mass, and the final BH masses, 
and the mass function of MBHs. In Section [4j we investigate 
the effects of Population III stars and the nature of metal 


poor globular clusters. Finally, in Section [5] we summarize 
our main conclusions. 


2 MODEL 


The starting point of our model is the formation of compact 
star clusters in major mergers (M 1 /M 2 < 3.5, Mi > M 2 ) 
of gas-rich disc galaxies. It has been shown in a series of 
numerical simulations that during such events gas initially 
in a rotationally supported disc looses angular momentum 
and collap ses toward the potential m i nimum of the merger 


remnant (Barnes & Hernquist 19961; Spring 

el et all 20051; 

Di Matteo et al. 20051: Naab et al. 

20061: 

Honkins et al. 

200(4 Cox et all 200(4 Li et al.1 2007 

: Maver et al. 201(tl. 


This gas will reach high densities on short time scales and 
be available to fuel star formation in the centre. The low 
amount of angular mo mentum facilitates the fo rmation of 
compact star clusters dRegan fc Haehneltl 120091 ). In prac¬ 
tice, the enhancement of mass inflow and the associated 
starburst due to the merger depend on physical conditions 
such as e.g. the mass rat io, impact parameter and inclina¬ 
tion dHopkins et al.ll2009f ). In this work, we consider major 
mergers with a parameter range leading to most efficient 
mass inflow and associated s tarburst, i.e., Mi/M 2 < 2 and 
an impact parameter < R v u (jKhochfar fc Burkertll2006 f. 

The progenitor disc is assumed to re late to t he host- 
ing dark matter halo via R d ~ XR v i r (e.g., iMo et aklll998l l. 
where A and R v i r are spin parameter and virial radius, re¬ 
spectively. In this work, we assume A = 0.05, however, our 
results are not depending on the specific choice of A as we 
will show below. 

We model the progenitor discs of merging galaxies as a 
Mestel, isothermal profile, i.e., £(r) = £o(-Rd/r), where £0 
and Rd are scal e parameters and estimated as a func tion of 
halo properties (|Mo et alJll99Sl : IDevecchi et alJl201(t) . 


£0 = 70 


14 \ 

15 km s -1 ) 


/disc \ 

005/ 



M 0 pc 2 , 


R d = 100 



R v 


700 pc 


pc. 


( 1 ) 


Here Vh is the circular velocity, and fdisc is the fraction of 
the disc mass with respect to the total mass including bulge 
and dark matter. In this wor k, we assume a fiduc ial value of 
A = 0.05 and /disc = A fe.g.. Ivan den Boschll200ll 'l. 

The amount of inflowing g as during the merger is esti¬ 
mated using the prescription in lHopkins et al.l (120091 ). which 
assumes that during the merger a gaseous and stellar bar de¬ 
velop which are out of sync and exert a torque on each other. 
The resulting mass inflow of gas is then estimated by, 


Minf = 27r£oRdI4(l — /gas )/disc'I'bar At, (2) 


where / gas is the gas fraction in the disc, 'Fbar is the mass 
fraction of stars in the bar, and At is the time since the 
merger. We assume /g as = 0.9, i.e., 10% of the disc is in stars 
before the merger takes place. The mass of the resulting star 
cluster increases with decreasing / gas , leading to formation 
of more massiv e BHs as we will show below. As suggested in 
iHopkins et alj (120091 '). we set \kbar = 1 based on our choice 
of mass ratio, inclination and impact parameter of mergers. 

We consider the initial 3 Myr after the merger, which 
corresponds to the life time of massive stars and the onset 
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Formation of Massive seed black holes 3 


of supernovae feedback which will halt star formation in the 
star cluster. 

Once mass inflow occurs the inner structure of the disc 
will change. However, we here focus on regions at galactic 
centres that are much smaller than these scales and we as¬ 
sume that the gas distribution is spherical. Prior to star 
formation we approximate the density profile of the central 
clouds with a singular isothermal profile, 

Cs 

Pgas = 2^ (3) 
where c s is the sound speed which we approximate using 
the virial temperature of the hosting halo. By integrating 
the density profile, we determine the radius of star clusters, 
i.e., M ini = ff cI 4nr 2 p gas dr, hence, 


r c i = 2.15 pc 


/ Cs \ 2 f Mi n { \ 

V 10 km s -1 / \ 10 5 Mg ) 


(4) 


In general, the star formation rate in such gas cloud can 
be parameterised by M star ~ V ^ gas where p < 0.05 is 
the star formation efficiency (e.g., iKrumholz et al.112012! ). 
The dynamical time of the gas clouds in our model is typ¬ 
ically < 10 4 yr which is much shorter than 3 Myr. There¬ 
fore, we assume the fiducial case of M c \ ~ Mint. On the 
other hand, before supernovae occurs, radiative feedback 
may reduce the conversion efficiency from gas to stars due 
to photo-evaporation of gas. Hence, we also study the case 
Of Mcl = 0.3Mi n f. 

The star clusters undergo core-collapse, and mas¬ 
sive stars migrate toward the centre of the star clus¬ 
ter. Recent N-body simulations show that the time scale 
of core-collapse is similar to that of dynamical friction 
dFuiii fe Portegies Zwartll20l4l . 


tcc ~ tdf 


1-91 r 2 cl a 

In (A) Gm max 


(5) 


where ln(A) is the Coulomb logarithm, a is the three- 
dimensional velocity dispersion of stars and m ma x is the 
maximum mass of stars based on the initial stellar mass 
spectrum. We assume that the stellar cluster is that of a 
King profile and use a Salpeter initial mass function with 
the mass range of m — 0.1 — 100 Mq. Due to core-collapse, 
the stellar density at the centre quickly rises. As a result, 
at the centre, stars frequently collide with each other and 
coalesce, resulting in t he formation of very massive star s 
(VMSs) of ~ 1000 Mq ( IPortegies Zwart fe McMillanll20C)3 ). 
Nu merical simulations estimate the fin al mass of such VMS 
as (IPortegies Zwart fe McMillanll2002h : 


My ms = m max + 4 x 10 -3 M c i/ c ln(A)ln 


f 3 Myr \ 

V to c ) 


( 6 ) 


where M c \ is the mass of the star cluster, and f c is the 
factor used to calibrate the analyt ical estimate against direct 
V-bo dy simulations (/ c = 0.2; IPortegies Zwart fc McMillanl 

mi). 


If the core-collapse time scale tec is longer than 3 Myr, 
stars will explode as type-II supernovae and we s et My ms = 
m max . Following IPorte gies Zwart fc McMillanl (120021 ) and 
iFuiii fc Portegies Zwartl ( 20141 ). we use a King profile with 
Wo = 3 for the stellar density distribution of our model star 
clusters. If the mass of the VMS s are greater than ~ 260 Mq, 
VMSs directly collapse to BHs dHeger fc Wooslevll2002l ) . Af¬ 
ter core-collapse, star clusters can shrink by a factor of a 


few dFuiii fc Portegies Zwartl 120141 ") , massive stars tend to 
distribute near the centres of star clusters due to mass seg¬ 
regation. We here will use the initial state of the star cluster 
in our calculations, noting that the calculated growth rates 
will be lower limits based on the fact that we do not take 
into account core collapse and mass segregation in the clus¬ 
ter. Subsequent disruptions of stars by the BH deplete the 
loss cone and drive the growth of the BH to a halt. How¬ 
ever, the loss cone can be refilled by stars that loose angular 
momentum and migrate towards the centre or have a high 
ellipticities. In this work, we consider angular momentum 
transport due to relaxation processes between stars. We will 
focus on two main types of relaxation processes, resonant 
and non-resonant. 

Two-body relaxation (non -resonant :NR) takes place 
over a time scales of (e.g., iBinnev fc Tremaind 120081 : 
iKocsis fc Tremai~nell201ll ). 


ttm(r) = 0.34 


G 2 pm2ln(A) ’ 


(7) 


where p is the stellar density, and m 2 =< m 2 > / < m > 
is the effective mass. For a Salpeter IMF int the mass range 
0.1 — 100 Mq, m 2 = 4.8 Mq. 

The second relaxation process we consider is resonant 
relaxation (RR). Stars close to MBH move on Kepler or¬ 
bits. These stars cause wire-like fluctuation in the gravita¬ 
tional potential governed by the MBH, and induce pertur¬ 
bations to stellar orbits. In particular, stars on the same 
plane exchange angular momentum via scaler RR resulting 
in a change of ellipticities of the orbits. As a result, stars 
with high ellipticity enter within the tidal radius of the cen¬ 
tral MBH and will be disrupted by it. The effect of scaler 
RR can be suppressed by Newtonian and general relativistic 
precession dHopman fc Alexander! f2006 ). in contrast to the 
classical Keplarian case. 

The scalar RR . tim e scale is estimated by 
( Hopman fc AlexandeiUhoQfjl l. 


= 


Arr 


N{< r) \ m 


(^] 2 p\r) 


1 

tm 


1 

tG R 


(8) 


wher e Arr is a numerical factor of order unity (~ 
3.56; iRauch fc Tremainel Il996h . S is a mean stellar mass, 
P(r) — 27r[r 3 /(GMBH)] 1 ^ 2 , £m and £qr are the Newto- 
nianand general relativist i c prec ession time scales. Follow¬ 
ing |Hopman_^Alexandeij (I2006T) . we estimate the Newto¬ 
nian and general relativistic precessions as follows, £m = 

A Ml ^^P(r),and£GR = f ( 4G m bh ) ^(r), where Am = 1 
is assumed. The shorter of the two relaxation time scales is 
used to estimate the loss cone refilling and growth of BHs, 


f relax (r ) = min(f N R, Irr) ■ 


(9) 


At radii where the interior mass is dominated by the 
stars and not the BH we only consider £nr, because the 
stellar orbits are not Keplerian anymore. Over the relax¬ 
ation time scale T^iax, the angular momentum of stars is 
transported, and hence, the loss-cone is refilled, causing the 
growth of BHs. We estimate the BH mass by integrating the 
stellar density pro file and time evolution as follows (see also 
IChen fc Liull2013h . 

M bh (£) = dt r Cl 47 rr2pstar(r ’ f) dr. (10) 

Jo Jo trelax \Tj t) 
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Since the accretion time-scale is shorter than the re - 
laxation time-scale of the star cluster (ISiokawa et alj|2015l) . 
we assume that stars entering the tidal radius are instanta¬ 
neously swallowed by BHs. Note that, however, this growth 
rate is an upper limit. In practice, some fraction of gas 
in stars can be unbounded when they are tidal disrupted 
dGuillochon fe Ramirez-Ruid l2013ll . Even initially bound 
gas, that settles in an accretio n disc, can be subject to feed - 
back from the accreting BH (IStrubbe fc Quataert ll2009ll . 
In addition, some stars near the central high-density re¬ 
gion can form binaries. If the binaries get close to the BHs, 
some stars can be kicked out due to three-body interaction. 
Thus, the growth rate may be reduced by factor 2-4, e.g. 
in cases where the mass fraction of the unbounded gas to 
the i nitial stell a r mass is ~ 0.5 and m out fl ow < m ac cretion 
dStrubbe fe Quataert II 2 OO 9 I '). 


3 RESULTS 

3.1 Formation of very massive stars 

Figure[l]shows the predicted mass of VMSs based on Eq.[6]as 
a function of the halo mass of the remnant. Most notably the 
mass of VMSs drops beyond halo masses of Mh ~ 10 8,5 ~ 9 . 
This is a direct consequence of the increasing core collapse 
time scale tec in larger haloes. The difference in predicted 
masses of VMSs as a function of redshift is again related to 
the longer core-collapse time scales at low redshifts, which 
result in systematic lower masses. Once the core collapse 
time scale is larger than the life time of massive stars, i.e. 3 
Myr, we assume that collisions between stars become negli¬ 
gible during the life-time of the star, and no VMS is formed. 
The long relaxation times would even in the case that a 
remnant black hole forms, be too long to grow it efficiently. 
The mass of VMSs increases with halo mass, because the 
mass of star clusters increases with halo mass in our model. 
However, with increasing cluster mass and radius, the core¬ 
collapse time becomes longer too. The mass growth of VMSs 
is logarithmic and only weak until tec ~ 3 Myr. Haloes with 
Mh ~ 10 9 Mq at 2 > 20 can host VMSs of ~ 1000 Mq. 
Note that, in this work, we assume low stellar metallici- 
ties (Z < 10~ 3 Zq). If stellar metallicities are higher than 
Z > 10 -3 Zq, stars experience strong stellar winds and loose 
a large fraction of their mass before their death. The max¬ 
imum mass of VMSs increases with redshift from 290 Mq 
for a halo of 1.1 x 10 s Mq at 2 = 7 to 1200 Mq for a halo of 
6.7 x 10 8 Mq at 2 = 20. In our model, M c \ for a specific halo 
mass does not depend on redshift, because in Equation [21 
Eo oc (1 + 2 ) 1 / 2 , Vh oc (1 + 2 ) 1 / 2 , and R v n oc (1 + 2 )" 1 . On 
the other hand, R c \ oc (1 + z) -1 , hence star clusters become 
bigger with decreasing redshift. 

As a result, the core-collapse time becomes longer, and 
the mass of VMSs decreases. If the stellar metallicity is 
l ow, even ~ 100 M q stars can directly collapse into BHs 
(lllegei- et al.l 120031 1. Feedback by supernovae can suppress 
further gas accretion or evacuate gas from galaxies. 

Once gas is evacuated from the disc and becomes part 
of the hot diffuse halo, it will take a cooling time, for it to 
join the disc again. Here we roughly estimate the feedback 
necessary to unbind the gas from the disc. The feedback en¬ 
ergy by SNe is E S n ~ 0.007 x 10 51 erg, where 0.007 



log M h (M e ) 

Figure 1. Mass of very massive stars in star clusters after 
core-collapse. Different colour lines represent different redshifts 
at which galaxies merge. Dashed line shows the threshold mass 
of stars that resu lt in black holes at end of their life time 
(iHeger et al.ll2003h . 


is the conversion factor from the star cluster’s mass to num¬ 
ber of massive stars causing type-II SNe given a Salpeter 
IMF in the mass range 0.1 — 100 M©. From Equation ([2]). 

QM 2 m 2 

-Esn oc Mh . The disc binding energy is E grav ~ R d oc . 
Then, EsN/E grav oc R V1I oc (1 + 2 ) -1 . Hence gas in discs can 
easily be evacuated at lower redshift, due to its lower den¬ 
sity. The threshold mass for haloes with E grav > Esn is 
M h > 1.6 x 10 9 (t^) 1 ' 5 - 

As shown in the figure 1, haloes of Mh ~ 10 7 ' 5 —10 9 Mq, 
which are in the typical mass range of the first galaxies 

dWise et al.ll2012bl : Ijohnson et alJl2013l : IPaardekooper et all 

l2013tl . form compact star clusters and cause the formation 
of VMSs. The VMSs of haloes with Mh < 10 7 ' 5 Mq result 
in pair-instability SNe dHeger et al]|2003h . 

3.2 Growth of BHs 

Once a VMS of > 260 Mq is formed collapse to a BH will 
ensue in a star cluster. In such a system RR/NR efficiently 
changes the angular momentum of stars, and refills the loss 
cone around the BH allowing it to feed. Figure [2] shows the 
RR and NR time scale as a function of radial distance from 
the centre of a star clusters. Due to general relativistic pre¬ 
cession a shortest time scale exists at a specific position. For 
a star cluster with M c \ = 10 5 Mq and R c i = 0.25 pc which 
is representative of the result of a merger between proto¬ 
galaxies with Mh = 3.8 x 10 8 Mq at 2 = 15, the RR time 
scale becomes shorter than NR at r < 0.02 pc. 

Hence, RR efficiently works for stellar accretion near 
the BH, and NR leads to the angular momentum change of 
stars in the outer parts. 
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As the stellar density decreases with time, the relax¬ 
ation time scale becomes longer. The RR and NR time scales 
at t = 1 Myr are shown as the dash lines in the figure. The 
time scales get longer by a factor > 2 from the initial state. 
The radius at which the RR time scale becomes shorter than 
the NR propagates outwards with time. In the fiducial case 
shown in figure [2] the radius is 0.1 pc at t = 1 Myr. Finally, 
depending on the RR and NR processes, the time scale for 
refilling of the loss cone can be shorter than ~ 10 s yr in the 
clusters. 

Following Equation 1101 the stellar density at a specific 
radius decreases with time due to RR and NR processes. 
Figure [3] shows the time evolution of the stellar density pro¬ 
file. Stars at r ~ 10~ 3 pc have short RR relaxation time 
scales and are quickly consumed by the central BH, which 
in turn generates a density gap. Stars close to the gap will be 
accreted next onto the BH. Initially the radius at which the 
total enclosed stellar mass is smaller than the central BH, is 
limited to r < 10~ 3 pc at t ~ 10 7 yr, therefore, all stars at 
large radii are only affected by NR. As the BH grows, the 
radius at which RR dominates becomes larger and reaches 
the radius of the cluster at t ~ 10 s yr. Note that, however, 
we do not follow the change of radial position of stars in 
this work. The actual stellar distrib ution can become mor e 
smooth in detailed simulations fe.g.. iFreitag &: Bendl200^1 . 
in which case the growth rate of the BH can be different. 

We follow the time evolution of the BH growth from 
the merging event till redshifts z = 6. The distribution of 
BH masses is shown in Figure [4] As shown in Figure [l] only 
haloes in the limited mass range {Mh ~ 10 7 — 10 9 M©) 
contribute to the formation of MBH seeds. 

As equation [7] shows, the NR time scale is proportional 
to the stellar velocity dispersion and the inverse of stellar 
mass density. Since the stellar density is almost constant 
despite of the cluster’s mass in our model, the BH growth 
rate is high even in low mass haloes or stellar clusters. How¬ 
ever, as it will be show below, most of the stars even in 
massive star clusters accrete onto a central BH by z = 6. 
Hence, the BH mass linearly increases with halo mass. We 
also compare the case without NR for galaxies that merge 
at z = 20, shown as the magenta dash line in the figure. 
Unlike in the case of both RR and NR processes active, the 
BH mass without NR is almost independent of halo mass, 
and smaller than that with NR. 

The radial dependence of the RR time scale implies that 
at small radii fast accretion of material onto the black hole 
will occur stalling the growth of the black hole until further 
material from larger radii is transported inwards. In the ab¬ 
sence of NR the transport of material from larger radii is 
very inefficient in effect limiting the black hole mass in our 
model to ~ 10 4 M©. Thus NR is able to tap into the mass 
budget of stellar clusters at large radii and RR efficiently 
funnels this material close to the black hole for accretion. 

Figure [3] shows the conversion fraction of the initial stel¬ 
lar mass of a cluster into a MBH by z = 6. As shown in 
Figure [2] the relaxation time scale is shorter than the aver¬ 
age duration from the first merging event to z = 6, hence 
most of the stars (> ninety per cent) accrete onto a central 
BH. With decreasing stellar density, the NR time scale gets 
longer, resulting in slow BH mass growth. The position at 
which the NR time scale becomes shortest gets closer to R c i 
as the cluster mass decreases, resulting in smaller conversion 



Figure 2. Time scales of stellar relaxation processes as a function 
of radial distance from centre of a star cluster. NR and RR rep¬ 
resent non-resonant relaxation (2-body) and resonant relaxation. 
Solid and dash lines show the relaxation time scale at t = 0 and 1 
Myr, respectively. The time scales are calculated for a star cluster 
with the mass of 10 5 Mg and the radius of 0.25 pc. 


factors of stars in clusters with > 10 5 Mg. In addition, as 
the stellar velocity dispersion increases with the star clus¬ 
ter’s mass, the NR time scale becomes longer, resulting in a 
somewhat inefficient conversion factor. 

In order to study the efficiency of stellar mass accre¬ 
tion, we compare it with the Eddington rate. Figure [b] rep¬ 
resents the time averaged mass accretion rate normalized to 
the Eddington rate. The mass accretion rate increases with 
redshift, because star clusters are more compact and the re¬ 
laxation time scale becomes shorter for mergers at higher 
redshift in our model. 

Due to the short relaxation timescale, the clusters lose 
most of their stars before z = 6. After most of the stars 
are swallowed, the stellar accretion rate rapidly decreases as 
shown in Figure 0 Therefore, the Eddington ratio at z ^ 15 
is lower than that at z ^ 10. 

However, even for star clusters formed at z = 25, the 
time average accretion rate is close to the Eddington rate 
(M/MEdd ~ 0.5). In addition, at Mh ~ 10 9 M©, the ac¬ 
cretion increases steeply. As shown in Figure [T| the initial 
BH (VMS) mass sharply drops at Mh ~ 10 9 M© due to the 
longer core-collapse timescale. On the other hand, the NR 
timescale is independent of the central BH mass. Hence, the 
Eddington ratio in the initial phase is very high. 

Thus, we suggest that the stellar accretion can be the 
dominant mode for the growth of BHs up to a mass of ~ 
10 5 Mg. For further growth of BHs to ~ 10 9 M©, efficient 
gas accretion is required. 
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log r/r 0 

Figure 3. Radial stellar density profiles. Different colours show 
the stellar density at different evolution times. The time scales 
are calculated for a star cluster with the mass of 10 5 Mq and 
radius of 0.25 pc. 



log M h (M 0 ) 

Figure 4. Black hole mass as a function of galaxy mass at z = 6. 
Different colour lines represent different redshifts at which galax¬ 
ies merge. The magenta dashed-line is the black hole mass without 
2-body relaxation for clusters formed at z = 20. 



Figure 5. Final BH mass at z = 6 normalized by the initial star 
cluster’s mass. Different colour lines represent different formation 
redshifts of the star clusters. 


3.3 Number density of seeds of massive BHs 

To estimate the cosmological relevance of the above dis¬ 
cussed formation channel, we combine our analytic model 
with a Monte-Carlo technique based dark matter merger 
tree. In combination with the dark matter mass function 
at any given redshift this will allow to predict the merger 
rate of galaxies and the associated black hole mass function 
(BHMF). 

Merger trees are constructed and app lied based on 
the extended Press-Schechter form alism (lLacev fe Cole 
1993J) using the algorithm presented in [Somerville fe Kolatt 
(1999). The merger trees have been su ccessfu lly imple¬ 
mented and tested in previous work (e.g.. iKhochfar fe Silkl 
120091 . 1201 lh . Our sample consists of 1000 merger tree reali¬ 
sations for each individual dark halo mass bin( 10 8 Mq to 
10 13 Mq.) Each tree is resolved down to a minimum halo 
mass of 10 7 Mq. We use cosmological parameters based 
on WMAP9 (Hinshaw et al. 2012), 12 a = 0.7, 12m = 0.3, 
as = 0.9 and h = 0.7. 

Figure [7] shows the redshift evolution of six sampled 
haloes which reach the halo mass of 10 s , 10 9 ,10 10 ,10 11 ,10 12 
and 10 13 Mq at z = 6 respectively. The yellow shade region 
represents the mass range in which the merging haloes can 
form MBH seeds. As seen in the Figure and expected, mas¬ 
sive haloes experience mergers in the shaded region at higher 
redshifts. For example, most haloes with Mh ^ 10 12 Mq at 
z = 6 can form MBH seeds at z > 15, resulting in plausible 
seeds of SMBHs at z > 6. In addition, due to the variation 
in the redhsift when the actual merger occurs, haloes with 
the same mass at z = 6 can host MBHs of different mass. 
Note that we assume the metallicity of the merging galaxies 
is low and remains so during the merger, and that the mass 


© 2016 RAS, MNRAS 000.H1HT1 




































Formation of Massive seed black holes 7 



Figure 6. Time-averaged stellar accretion rate normalized by the 
Eddington rate. Different color lines represent different redshifts 
at when galaxies merge. 



10 15 20 

redshift 


Figure 7. Example redshift evolution of haloes. The yellow 
shaded region indicates the mass range in which haloes can form 
BH seeds when they experience major merging. 


loss of VMSs is not significant. If galaxies are metal-enriched 
to a high degree, the initial BH mass is likely to change. 

Major mergers with mass ratio <1:3 are effective 
means in transporting angular momentum and allowing the 
gas to collapse to the central region of a galaxy (e.g. Bois 
et al. 2010). The efficiency of this process decreases with 
increasing mass ratio and we will in the following focus only 
at major mergers with mass ratio <1:2. We have checked 
that inclusion of minor mergers does not change the BHMF 
significantly. We identify the first major merger event in the 
history of a halo as the BH seeding event. 

Figure [8] shows the predicted BHMFs at z = 6,10,15 
and 20. At 2 = 20, the BHMF is limited to Mbh < 10 4 Mq 
because the merging events just took place. As time pro¬ 
gresses the number density of BHs and their average mass 
increases. 

At z ~ 10 — 15 we find MBHs of ~ 10 5 Mq with 
number densities of ~ 10~ 6 to 10 - ' Mpc~ 3 . If we assume 
M c \ = 0.3M inf , the number densities decrease to ~ 10~ 8 
to 10 -9 Mpc -3 which are simila r to that of observed QSOs 
at z > 6 (e.g., iFan et ad] l2006l l . This is because the core¬ 
collapse time scale gets longer, hence the halo mass range 
becomes narrower. In addition, the growth rate of MBHs by 
stellar relaxation also becomes lower. 

MBHs of 10 s Mq that grow close to the Eddington limit 
will reach ~ 10 9 Mq by z ~ 6 — 7 and thus could be the 
potential progenitors. To illustrate this point, one can con¬ 
sider a MBH formed at z = 15 (10), it can reach 10 9 Mq at 
2 = 6 with an average Eddington rate of 0.6 (0.9). 

Cosmological simulations show that 10 s Mq MBHs at 
galactic centres can grow to ~ 10 9 Mq by z ~ 6 via 
gas accretion using a the Bondi-Hoyle model at the sub - 
grid level (e.g, iDi Matteo etU] 120051 . 120121 : [lT et al. !ioo3). 


These simulations show that once a MBH is formed at the 
centre of galaxies in the halos with Mh > 10 10 Mq, the 
gas keeps accreting onto the MBH close to the Edding¬ 
ton li mit even un der the pr esence of feedback (see however, 
Ijohnson fc Khochfarl [201ll . for high resolution simulations 
including radiative transfer.). If we allow the MBH seeds of 
10 5 Mq to grow at the Eddington limit, the time needed to 
grow from 10 s to 10 9 Mq is t — yo^TsaihilO 4 ~ 0.46 Gyr 
for e = 0.1, where rs a i is the Salpeter time scale = 

Hence, in order to explain SMBHs of 10 9 Mq at 2 ~ 6 (7), 
MBH seeds of 10 5 Mq are required at 2 > 10 (13). Our 
models satisfy these constraint. 

It is well known observationally that the mass of MBHs 
at galactic centres in local galaxies ti ghtly correlates with the 
stellar mass of galac tic bulges fe.g.. IMarconi fe Huntll2003l : 
IGiiltekin et alj|2009f l. On the other hand, at high redshifts, 
the BH-stellar mass relation is less well constraint, due to 
difficulties of observing BHs of < 10 s Mq. If we compare our 
modeled BH mass density with recent obser vational results 
of stellar mass densities rtStark et alj|20l3 l. the BH mass 
density is approximately 5 orders of magnitude smaller than 
the observed stellar mass density at 2 ~ 10, and continues to 
become even smaller with respect to the stellar mass density. 
This implies that gas accretion potentially plays a dominant 
role in the growth of BHs after MBH seeds of ~ 10 5 M 0 
form. 


4 DISCUSSION 

4.1 Population III stars 

In this work, we consider a Salpeter IMF in the mass 
range 0.1 — 100 Mq. However, the metallicity of gas is 
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Figure 8. Black hole mass function at z = 20,15,10 and 6. Solid 
and dash lines represent the case of esF = 1 (our fiducial model) 
and 0.3, where egp is the conversion efficiency from gas to stars 
in the inner region r < r c \. 


low at high redshift, and hence the IMF can tend to be 
top-heavy. As an upper limit, we consider zero metalic- 
ity stars, i.e., Population III (Pop III) stars. The IMF 
of Pop III stars is still unknown due to missing obser¬ 
vational probes as well as uncertainties in the modelling 
of their formation. Given that accretion rates onto proto- 
stellar seeds in primoridal gas clouds are high it is gener- 


their metal enriched counter parts (e.g., Omukai & Nishi 

19981; Omukai & Pallal 120031; iNakamura & Umemurd 

2001 ; 

Bromm et al. 20021: Abel et al. 2002: Yoshida et al. 

2006, 

2008 

Turk et al. 12009 

: Stacv et al.. l20ld. 2012- Clark et al. 

2011 ; 

Umemura et al.l 

20121: Hirano et al.l [2014t: ISusal 20131: 

Susa et al. 20141). We 

re-examine our model trying to ac- 


count for such massive Po p III stars assumin g a top heavy 
IMF from 10 — 500 M@ IlHirano et al.1120141 ) and a slope 
a = —2.35. Pop III stars usually for m in mini-haloes 
of ~ 10 6 Mfl (e.g., iBromm et al.l 20021). as a single star 
(lYoshida et al]|200a) or binary I Turk et al.ll200^ ) or multi¬ 
ple ( Clark et all 201 lh . On the other hand, UV radiation in 
the Lyman-Werner band from star forming g alaxies can sup¬ 
press star formation in neighbouring haloes (I Johnson et al.l 
120131) . As a result, there may be the possibili ty that atomic 
cooli ng haloes host massive pristine gas (e.g.. lAgarwal et al] 
Ell). The resulting BHMF at a = 6 is shown in Figure O 
In contrast to star clusters with a fiducial IMF star clusters 
with a top-heavy IMF produce more massive black holes at 
their centre per unit mass in stars. Massive Pop III stars 
cause core-collapse of more massive haloes within 3 Myr, 
and lead to the formation of heavier VMSs. In addition, be¬ 
cause of the increase of the average stellar mass < m > from 
0.3 Mq for the normal IMF to 30 Mq for the Pop III IMF, 
the relaxation time scale becomes shorter by factor ~ 10 


(Equation [7] and [8j. As a result, massive BHs of ~ 10 6 Mq 
can form by 2 ~ 6. 

In other words, for such a top heavy IMF, even star 
clusters with similar size and density of typical globular clus¬ 
ters in the local universe can cause core-collapse, and form 
VMSs. 

Note that, however, that massive galaxies (> 10 s Mq) 
are li kely metal enriched due to type-II SNe of prior stars 
(e.g.. lMaio et al.1 l201ll ; IWise et al.ll2012al ). Therefore, mas¬ 
sive star clusters consisting of Pop III stars may not form in 
practice. On the other hand, if metal mixing is not efficient, 
some pockets of pristine gas will be still able to lead to the 
formation of Pop III sta r clusters, with sh ort core-collapse 
time scales (see however. iSmith et al.ll2015l ). 

Local metal poor globular clusters can give constraints 
on the stellar IMF, density and the cluster’s size to explain 
their old-age of > 10 Gyr. Based on our results star clusters 
with a local IMF and initial size of > 1 pc should be long- 
lived. 


4.2 Metal poor globular clusters 

Merging haloes between Mh ~ 10° — 10 10 Mq in our model 
result in the formation of globular clusters. We here present 
the mass function focusing on those with tec > 3 Myr and 
J?sn > -E/gra V ■ For the star clusters with tec > 3 Myr, central 
massive BH seeds via core-collapse do not form, leaving the 
stars unaffected in the cluster. In addition, if Esn > -E grav , 
SN feedback is likely to suppress further gas accretion onto 
the star clusters, resulting in isolated systems that are dense 
and compact. It is widely known that observed globular clus¬ 
ters do not have dark matter haloes. On the other hand, 
some theoretical models showed the formation of globu¬ 
lar clusters in dark matter haloes, and the number density 
nicely fit with observatio ns dMoore et al .IIJOOgI ; iTrenti et al.l 
l2015tlKimm et ahll2015l) . These authors assumed the host¬ 
ing dark matter haloes would subsequently be stripped away 
by tidal interaction with galaxies. Although the formation 
mechanism of globular cluster is still unclear, we here present 
the predicted mass distribution from our merger model as a 
possible sub-dominant channel to forming globular clusters. 

Figure QT)] shows the mass function at 2 = 6. The mass 
function has a peak at M ~ 1.5 x 10 5 Mq and the mass range 
is consistent with that of observed globular clusters in local 
Universe. Recent observation of metal-poor globular clus¬ 
ters (MPGCs) show that the stellar popu lation is very o ld, 
and c ompatible with formation at 2 > 6 (IBrodie fe Straden 
120061) We compare the number density of MPGCs in our 
model with that of the observed one in local Universe. Thus, 
our model alone cannot reproduce the number density of the 
observed local MPGCs, while it matches their physical prop¬ 
erties (size and mass). 

The subsequent interaction with their envir onment via 
proce sses like tidal stripping and merging (e.g, iLeon et al.l 
l2000l) will change the mass function of MPGCs and its final 
outcome will depend on th e futur e evolution. 

Recently ITrenti et al.l (120151) modeled local MPGCs 
in co smological N-body simulations (see also iMoore et al.l 
120061) . They assumed GCs form in merging primordial haloes 
of 10 8 Mq at 2 ~ 10, and used their model to explain the 
observed nature of local MPGCs, e.g., age and metallicity 
distribution, and their spacial distribution. Our model of GC 
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Figure 9. Black hole mass function at z = 6 assuming a top 
heavy IMF in the mass ranges 10 to 500 Mq with slope a = 
—2.35. Dash lines show our fiducial model (M = 0.1 — 100 Mg) 
as shown in Figure [8] 


formation requires somewhat more massive merging haloes 
(> 10 9 Mq) suggesting a possible halo mass bias leading to 
the formation of black hole seeds. 


5 SUMMARY 

In this work, we analytically model the formation of MBHs 
in the first merging galaxies. These BHs provide the natu¬ 
ral seeds for SMBH at z ~ 6. We show that compact star 
clusters that likely form during such mergers are prone to 
core-collapse and produce VMSs of ~ 1000 Mq at their cen¬ 
tres. VMSs collapse to BHs without going through a SN 
stage and subsequently grow via swallowing stars in their 
vicinity. Relaxation processes are efficient in refilling the loss 
cone around the BH continuously with new stars. We find 
that the BHs mass sensitively depends on the the radius 
and mass of the hosting star clusters which determine the 
core-collapse and stellar relaxation time scales. Within this 
scenario major mergers of galaxies with > 4 x 10 s Mq at 
2 ~ 20 lead to the formation of > 10 5 Mq BHs by z ~ 10 
which are likely progenitors of SMBHs at z > 6. 

Based on our results the average relation between host 
stellar mass and black hole mass will not be the same 
as observed in the local universe. After the formation of 
first massive BHs, the BHs can grow by gas accretion, and 
their mass can be clo se to the locally observed relation 
(iDi Matteo et al1l2005l ). In addition, such deviation will also 
reduce over time as merging of galaxies and black holes will 
move host galaxies and black holes closer to the locally ob- 
served relat i on as a consequence of th e central limit theorem 
l|Pendl2007l ; iHirschmann et aI1l2010l) . 

Our model predictions are sensitive to the properties of 
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log M gc (M 0 ) 

Figure 10. Mass function of globular clusters resulting from 
mergers of galaxies. 


the hosting star clusters, e.g., mass, radius, IMF of stars. We 
here assumed physically motivated analytic relations that we 
will further test and systematically investigate in a follow-up 
study using detailed numerical simulations. 
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